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Abstract 

A spectroscopic method for staggered fermions based on thermodynamical considerations is proposed. 
The canonical partition functions corresponding to the different quark number sectors are expressed in 
the low temperature limit as polynomials of the eigenvalues of the reduced fermion matrix. Taking the 
Oh' zero temperature limit yields the masses of the lowest states. The method is successfully applied to the 

Goldstone pion and both dynamical and quenched results are presented showing good agreement with 
that of standard spectroscopy. Though in principle the method can be used to obtain the baryon and 
dibaryon masses, due to their high computational costs such calculations are practically unreachable. 
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One of the most critical steps in hadron spectroscopy is the choice of the wave function for the particle of 
interest. Without finding the proper interpolating operator contamination from other particle states can 



often occur. This does not cause serious difficulties when one deals with simple particles such as mesons 



o 
o 



or baryons, but the extraction of the correct mass can be a tremendous task when one is looking for more 
compound states. Therefore, a spectroscopic method that does not require knowledge of the wave function 
is desirable. 

The question was first addressed in Ref. [1], where the explicit inversion of the fermion matrix on dupli- 
cated lattice configurations and the examination of the exponential decay rate of the lowest hadron state led 
$_i ' to the expression 

am-,, = • max In | | 2 (1) 

Lt |A fe |<i 

for the mass of the Goldstone pion on each configuration, where a is the lattice spacing, L t is the number 
of lattice sites in the temporal direction and are the eigenvalues of the reduced staggered fermion matrix 
(see Section 3). The aim is to find the relation between the hadron spectrum and these eigenvalues. 

There has been many advances recently in the canonical approach to finite density QCD [2-5]. Based 
on the canonical formulation we make an attempt to clarify and extend the findings of Ref. [1] and give a 
method which in principle can be used to obtain the masses of different particles. 

The paper is organized as follows. In Section 2 we summarize how canonical partition functions can be 
used to obtain the masses of particles. This is followed by a description of how canonical partition functions 
can be obtained on the lattice in Section 3. The way one can find the relevant eigenvalues is shown in section 
4. In Section 5 we show how the Z 3 symmetry can be used to simplify our formulae. The case of baryons is 
explained in Section 6 while the case of mesons is discussed in Section 7. Finally, after showing our numerical 
results in Section 8 we conclude in Section 9. 
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2 Masses from canonical partition functions 

Let us consider the general case when we have n s different quark fields. Let Ni and \ii denote the quark 
number operator and the quark number chemical potential corresponding to the zth quark field, respectively. 
Then the grand canonical partition function at a given set of chemical potential values (/ii , /j,2, ■ ■ ■ , Hn a ) and 
temperature T is given by 



Z(mi,M2, • ■ • ,Vn„T) = Tr 



-(H-mNi-^^- ... -ij. ns N ns )/T 



(2) 



The canonical partition function corresponding to a given set of quark number values N\,..., N Ug can be 
obtained by taking the trace only over the subspace N± = N± , . . . , N ng = N ng . 

Z Nl ,...,N ne (T) = Tr 

= Tr 
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When one introduces imaginary chemical potentials [6], the different canonical partition functions become 
the coefficients in the Fourier expansion of the grand canonical partition function. 

i ,-2kT ,-2-nT 

Z Nl ,., Nn AT) = — J o d/ii...y d f i ns e- i ^ T ---e- i ^ N "°' T Z(m,...,i» ns ,T) (4) 



Z(i^,...,i^ na ,T)= ^ ... Y, Z Nl _ Nns (T)e i ^/ T ---e i ^nJT (g) 

N±=— oo N ng = — oo 

When the aim is to find the energy of the lowest state in the sector corresponding to quark numbers 
(Ni , . . . , N na ) one has to examine the low temperature behavior of the free energy 



F Nl ,..., Nna (T) = -TlnZ Nu ..., Nna {T). 
The canonical partition function can be written as 

oo 

7 m _\^„( Jv i.-. JV ».) P -Bi* ri --* r "' ) /T 

A/Vi,...,JV„ s (.-t ) - n k e 
fc=0 



(6) 



(7) 



respectively. In sector (0, . . . , 0) the lowest state is the vacuum, which is assumed to be non-degenerate. 



where ^ 1 ' -' " s; and 

^ivi,...,iv„ s ,i ^ e energy and the multiplicity of the kth state in sector (N\, . . . , N ns ), 



(8) 



k=i 



Then the difference of the free energies of sector {N\ , . . . , N ns ) and the zero quark number sector can be 
rewritten as 

F Nl ,..., Nna (T) - F ,...,o(T) = 2#*.--"».) - - Tln4 JVl '-' JV ' lJ 

00 (ATi,...,JV ns ) 
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The mass of the lowest state in sector (N\ , . . . , N ng ) is the difference of the energy of the ground state in 
this sector and the energy of the vacuum state, 



(JVi, 



E (»U...,»n.)_ E W,..,0)^ (1Q) 

If the temperature is much smaller than the energy differences that appear in the exponentials in equation 
(9) then the last term on the r.h.s. of (9) is negligible compared to the other terms. In this region the difference 
of the free energies follows a linear behaviour, 

fiV!,..,^,) 



F Nl ,..., Nn ,(T) - F 0i ... t0 (T) 



rlnn W,.,A^J 



(11) 



where the slope of the linear behaviour depends only on the multiplicity of the ground state. Therefore, the 
mass of the lightest particle carrying quantum numbers (Ni , . . . , N ns ) and its multiplicity can be obtained 
by a linear extrapolation to the T = limit. 



mi Nl '-> N ^ = Urn [F Nu ..., Nne (T) - F ,..., (T)] 



(12) 



3 Canonical partition functions on the lattice 

The temperature on the lattice is given by T = l/aL t , where L t is the number of sites in the temporal 
direction and a is the lattice spacing. Let /tj = denote the chemical potentials in lattice units. In order 
to introduce these chemical potentials on the lattice the forward time-like links have to be multiplied by e 1 ^ 
and the backward time-like links by e~ 1 ^ in the fermion determinant of quarks of type i [7]. Then the grand 
canonical partition function using staggered lattice fermions can be written as 



Z(»£i,... ,*£*.)= J[dU]e- s »M l[det Mim^t^Uy 



./4 



(13) 



where m, denotes the bare mass and rij denotes the number of tastes of the ith staggered quark field. The 
functional integral is taken over all possible gauge configurations U and S g [U] denotes the pure gauge part 
of the action. 

The partition function can be rewritten in the form 



Z(ip,i, ...,ip, 



»J = / [dU] 



e ]JdetM(mi,0,U) ni ^ x 



Sg[U] 
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det M(m,i, ifc, U) \ n ^ 4 
dct M(nn, 0, U) 



(14) 



»=i »=i 

The ratios of the determinants can be treated as observables while the functional integral can be taken using 
the measure at /tj = 0. Then the partition function becomes the expectation value of the determinant ratios 
taken over the ensemble generated at zero chemical potentials, 



Z(ifn, . . .,i/t„J = z ■ (Yi 



det M(m,, ifa, U) 
det M(m 4 , 0,[/) 



(15) 



where Z denotes the zero chemical potential value of the partition function [8] . 

Therefore, the canonical partition functions are obtained by taking the expectation values of the Fourier 
components of the determinant ratios. 
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In order to perform the assigned Fourier transformations we need the analytic /i-dependence of det M[ij]b). 
In temporal gauge, the fermion matrix can be written as 



M{ijl) 
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(17) 



B Lt -iJ 
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where U denotes the remaining time direction links on the last timeslice (including the correct staggered 
phases) and Bk is the spacelike staggered fermion matrix on timeslice k. In matrix (17) each block is a 
3V x 3V matrix, where V = L 3 and L s is the spatial size of the lattice. After performing L t — 2 steps of 
Gaussian elimination, the determinant of (17) can be written as 



det M{ifi) = e 3VLt ^ det (S - e"^*) , 



where 
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(18) 



(19) 



is the 6V x 6V sized reduced fermion matrix [9]. Let Afe denote the eigenvalues of S. Then (18) can be 
written as 



6V 



and thus, 



det M(ifi) = e 
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The eigenvalues of matrix 5 have a symmetry, according to which whenever A is an eigenvalue of S then 
1/A* is also an eigenvalue of S [9]. Therefore, each eigenvalue whose absolute value is greater than 1 has a 
pair with an absolute value smaller than 1, and vice versa. (We will not deal with the case when at least one 
of the eigenvalues lie on the unit circle because these gauge configurations constitute a zero measure set.) 
Then (21) can be written as 
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(22) 



where the product is taken over only the eigenvalues lying inside the unit circle. From now on when the 
limits of a sum or product taken over the eigenvalues of S are from 1 to 3>V , then the sum or product is 
meant to be taken over only the "small" eigenvalues, that is, the eigenvalues with absolute value smaller than 
1. 

When the temperature is low (T <C T c ) a gap appears between the "small" and "large" eigenvalues of S 
(see Figure 1). This makes a Taylor expansion of (22) in the small eigenvalues possible. As the temperature 
decreases the small eigenvalues become exponentially smaller, increasing the validity of the series expansion. 
Including the rational exponent for the number of tastes nt a first order expansion gives 
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(23) 



By performing an nth order Taylor expansion we explicitly obtain all the Fourier coefficients up to nth order. 

This way the assigned Fourier transformations in equation (16) can be performed easily configuration by 
configuration by simply choosing the coefficients of the corresponding exponential terms. The order of the 
leading order term for sector (iVi, . . . ,N ns ) is \Ni\ + ■ ■ ■ + \N ns \. When all the quark fields have 4 tastes 
(rii = 4) the leading order term for sector (Ni , . . . , N ng ) can be written as 
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(24) 



where *(sgn Ni) in the exponent means that there is a complex conjugation if N is negative. A fc stands for 
the kth eigenvalue of the reduced matrix obtained from the fermion matrix of the ith quark field. 
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Figure 1: The absolute values of the eigenvalues of the reduced fermion matrix 5ona typical 6 3 x 24 lattice 
configuration at a temperature of T « 25 MeV. At this temperature there is already a noticeable gap between 
the eigenvalues lying inside the unit circle and the ones lying outside. The eigenvalues that are relevant for 
calculating the canonical partition functions are circled. 



The leading order term in case of arbitrary number of tastes can be obtained from (24) as follows. The 
formula within the expectation value signs can be written as a homogeneous polynomial of the eigenvalues 
of degree |i\Ti| + ■ • • + \N n , | using the expressions 



^ (^fe^) j = 1 , . . . , iYj if Ni is positive, and 

fc=i 

3V ^ . 

( A fe°*) j = 1, . . . , -Ni if N is negative. 



(25) 



k=l 

The leading order term in the general case is obtained by replacing the expressions 

E(^y ^E(^y and £(^t with t£(^t ^ 

k=l k=l k=l k=l 

in the above polynomial. 



4 Obtaining the relevant eigenvalues 

In order to calculate the canonical partition functions using the description given in (24), (25) and (26) we do 
not need all of the 3V small eigenvalues. At lower temperatures (L t > 50 — 100) the small eigenvalues alone 
span a range of 20-40 orders of magnitude. Therefore, the relevant eigenvalues that contribute significantly to 
the sums in (24)-(26) are the largest few of the small eigenvalues (see Figure 1). Since the condition number 
of matrix S at low temperatures can be in the range of O(10 60 — 10 100 ) and these relevant eigenvalues are 
in the middle of the spectrum, finding these eigenvalues seems practically impossible. 

Nevertheless, the matrix S has some symmetry properties that make it possible. The spacelike staggered 
fermion matrices Bk, which appear in (17) and (19), obey a 75-hermiticity 

l5 B k =B[ lb , (27) 

where 

M xv = S xy -(-l)^ x ". (28) 
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Figure 2: The absolute values of the relevant eigenvalues found by the double precision version of ARPACK 
(red crosses) and their relative errors (green squares) on a typical lattice of size 6 3 x 100. 

Therefore, the inverse of S can be obtained as 

= <-i> w {I 7 S ) S, (° "?)• < 29 » 

As a consequence, once we have the matrix S both S + S^ 1 and S — 5 _1 can be easily constructed. Then 
by inverting these two one can arrive at 

Q= 1 -[(S + S-r 1 -(S-S-r 1 ]- (30) 

The order of magnitude of the condition number of S + S^ 1 and S — S^ 1 is less than half of that of S. 
Therefore, much less numerical precision is sufficient for their inversion. 

If Afe is an eigenvalue of S then Afc/(1 — A|) is an eigenvalue of Q. If Afe is a small eigenvalue, then 
|A|| <SC 1. In this case using Afc/(1 — A|) for the calculations instead of A& does not make any difference. If 
Afe is a large eigenvalue of S, then Afe/(1 — A|) « — 1/A|, which is negligible compared to the relevant small 
eigenvalues. That is, the relevant eigenvalues of S become the largest eigenvalues of Q. 

The procedure for finding the relevant eigenvalues was as follows. After fixing the temporal gauge the 
matrices Bk were built. Then from (19) and (29) the matrices S and S~ x were constructed. Since S and 
S^ 1 are very badly conditioned, their construction as well as working with them requires high precision. 
For these and the latter calculations the GNU multiple precision arithmetic library (GNU MP) was used. 
Then after inverting S + S^ 1 and S — S^ 1 the largest several eigenvalues of Q (the relevant ones) were 
obtained using the double precision version of ARPACK. This last step may sound doubious but in fact the 
double ARPACK was found to be able to reliably find the eigenvalues that were not more than 10 orders of 
magnitude smaller than the largest one (see Figure 2). 

5 Consequences of the Z 3 symmetry 
5.1 Consequences for ^jv 1 ,...,jv„ s 

From (22) it can be seen that the quantity within the expectation value signs in equation (15) is periodic 
in each fa with a periodicity of 2ir/L t configuration by configuration. Therefore, the lattice grand canonical 
partition function (13) is also periodic with 2-k jL t in each /V 

Performing a Z3 transformation, that is, multiplying all the time-like links on the last timeslice of an 
SU(3) configuration U by ej fa = e 2 ™ J/ 3 , j = 0, 1,2) gives another SU (3) configuration denoted by U ej . 
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Then the partition function can be written as 



1 2 n 3 

Z(iAi,...,t/i n .) = -y; /[dC^e"^ 'I []dctM(m 4 ,zA„^r /4 - 

6 3=0 7 i=l 



(31) 



The functional measure and the gauge action are both symmetric with respect to Z3 transformations [10]. 
Thus, 

Z(ijl 1 ,...,ifi ne ) = AdC/le-^^i^fjdetMK,^,,^)"* 74 

J j=0 i=l 

= / [dU]e- s °M ndctMK,0,t/r/ 4 x I Vft fdetM^ifc.C^V'' 4 



^A±V detM(m;,0,£/) 



(32) 
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7—0 2—1 v ' 



Since the {/ — > J7 £j transformation can be applied in eq. (18) by simply multiplying S by £j, the ratios 
of the determinants in (32) can be rewritten as 



detM(m,i,ifli,U £ i) 
detM(m i; 0, U) 

Combining (33) with (32) we obtain 
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(34) 



which means that the grand canonical partition function has an extra periodicity: if 2ni/3L t is added to all 
the chemical potentials then the value of the partition function remains unchanged [6]. 



. 2tt .„ . 2tt 

z{i/j ll ,...,ifj, ne ) = z\ m x +i—,...,inn e +l ^i t 



(35) 



Therefore, the canonical partition functions Zjv^...^ where the total number of quarks N\ + ■ ■ ■ + N na is 
not divisible by 3 are zero [4]. 

Taking this into account the expectation value of the first order expansion in (23) gives 



dot 



detM(O) J 



nt/i\ I 3V 3V 



(36) 



fc=l 



fc=l 



5.2 Application on a term by term basis 

Let A[U] be a gauge invariant quantity (a gauge invariant function of the gauge configuration U). Then the 
expectation value of A[U] is 



{A[U}) = ^ j [dI/]e- s «[ l/ ldetM(mi,0,^) ni / 4 ••• det M(m n „ 0, [/)™" s/4 • A[U]. 



(37) 
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Using the Z% invariance of the gauge action and the integration measure one can rewrite (37) as 



(A[U}) = ± J[dU], 



s a [u] -Q det m(TOjj o, £/)"* /4 x - 52 A[U e 



i=i 
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dctM(m t ,0,L^) 
det M(m,i,0,U) 



i/4 



that is, 



j=0 



i=l 



det M(mi, 0,17^) 
detM(m l7 0, £/) 



,/4\ 



(38) 



(39) 



Using (33) an expansion similar to (23) can be applied to these determinant ratios. 

All the quantities of the form of (25) are gauge independent, therefore, each term of the series expansion 
can individually be taken as A[U]. This way the Z3 symmetric form of (39) can be applied to each term in 
the series expansion. As an example, if we have only one staggered field (n s = 1) with number of tastes rit 
then the expectation value of the first order terms in eq. (23) up to leading order become 
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Applying this technique term by term the series expansion of (23) and (36) up to third order becomes 
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(41) 



In the third order expansion in (41) all the terms are already Z3 invariant. If a term in the series expansion 
is Z3 invariant then it does not change when the procedure of (39) is applied to it. When the procedure 
(39) is applied to a non-Z3 invariant term, its expectation value becomes the expectation value of the sum 
of higher order terms. This procedure can be continued order by order and as a result, all the remaining 
terms in all the quark number sectors of the series expansion of Z(ifii, . . . , ifi ns ) become Z3 invariant. 



6 Application to baryons 

In principle the method described in Sections 2 and 3 can be used to measure the mass of the lowest state 
in any quark number sector. For example, one can think of the di-baryon (deuteron), or the bound states 
of even more baryons. However, in most cases technical difficulties occur. Let us examine the case when 
one tries to measure the mass of a baryon, for example the proton. For that we use two staggered quark 
fields, one for the u quark with n u tastes and one for the d quark with rid tastes. (We omitted the third 
light quark, the s quark, the inclusion of which in our case does not change the picture significantly.) 

The proton is believed to be the lowest state in the N u = 2, Nd = 1 channel, therefore, according to (11) 
we need to examine the low temperature behaviour of 

F Nu=2 . Nd=l {T) - F N ^ Nd=0 (T) = -Tin (f^j) • (42) 
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According to (24)-(26) the relevant canonical partition functions to leading order are 

z - Sz -(^(l^(l^)-w(|^) 2 (|^)>^ m 

Zo,o = Z-(l), 
therefore, for the proton mass one obtains 

— A-H^ffi^ffi^-^fi^'S^))- 1441 

As the temperature decreases (L t increases) the eigenvalues become smaller and smaller, and only the leading 
order term matters in the limit. 

The formulae for the masses of the 2-baryon, 3-baryon, etc. channels can be obtained similarly. These 
can in principle be used to measure the bonding energy of several-baryon states. 

The description seems simple, but there is one difficulty. The expression of which the expectation value 
is taken in equation (44) can be any complex number, whose real part can be both positive and negative. 
Its expectation value is much smaller than its value at a typical gauge configuration. On a 6 3 x 24 staggered 
lattice with a s=s 0.33 fm, ro, w 330 MeV and T w 25 MeV this value at a typical gauge configuration is of 
O(10~ 10 ) while the expected order of magnitude of the expectation value is O(10~ 20 ). That means that the 
number of configurations needed for a correct result would be of O(10 20 ), which is prohibitive. The problem 
becomes even more severe when one decreases the temperature in order to get closer to the T — > limit. 



7 Application to mesons 

When nd = n u = n t /2, nid — m u and we are looking at one of the Nd — —N u sectors this sign problem 
does not arise. These sectors can be labelled with one parameter, the third component of the isospin 
h = {N u - N d )/2. Since = \f ] for all k, we will write X k only. 

The lowest state in the h = 1 sector is expected to be the Goldstone pion. Its partition function can be 
written as the expectation value 



lo „ 

= 1 = ^N u 



--l,N d = 
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(45) 



which is a manifestly positive polynomial of the eigenvalues. Therefore, it can be easily evaluated, and by 
taking the zero temperature limit 



ami 3 



1 

Lt 



64 



3V 

fc=i 



(46) 



one directly obtains the mass of the lowest state in the I3 — 1 channel. 

The formulae for the energies of the lowest state in higher I3 channels can be obtained similarly. These 
can be used to investigate pion-pion scattering and several-pion states. 

The result for the pion mass given in equation (46) obtained using purely thermodynamic considerations 
can be compared to formula (1). If L t is large, then after taking the logarithm the factor n 2 /64 gives a 
negligible contribution compared to that of the sum. If the temporal extension is large even compared to 
the spatial volume, then the sum in (46) is dominated by the largest of the small eigenvalues. In this case 
equations (46) and (1) evaluated on a single configuration yield approximately the same results. However, 
while Ref. [1] only states that relation (1) holds configuration by configuration and does not mention how to 
obtain results over an ensemble of configurations, equation (46) describes a method for taking the ensemble 
average. 
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50 


331 


322 




100 


1196 
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200 
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467 


300 
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255 





Table 1: The number of configurations used for dynamical staggered calculations with a spatial volume of 
6 3 . 
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aFj 3 —i — aFj 3 - 


=0 




n t = 2 


n t = 4 


n t = 8 


50 


0.5344(12) 


0.4971(12) 




100 


0.5066(2) 


0.4826(4) 


0.4639(4) 


200 


0.4931(2) 


0.4760(1) 


0.4641(3) 


300 


0.4876(3) 


0.4730(3) 




— > oo 


0.4787(3) 


0.4688(3) 


0.4643(7) 




0.47864(3) 


0.46903(4) 


0.46426(3) 



Table 2: The differences of the free energies, their T — > extrapolated values and the spectroscopic pion 
masses on dynamical staggered configurations with a spatial volume of 6 3 . 



8 Results 

8.1 Dynamical staggered fermions 

We performed calculations using dynamical staggered configurations to measure the masses in the first 
isospin channel as described in Section 7. We used the Wilson plaquette action for the gauge fields and 
unimproved staggered fermion action. In order to be able to check whether the root taking of the fermion 
determinant changes the results significantly, calculations were done using rooted staggered fermions with 
n t — 2 (n u = nd = 1) and nt = 4 (n u = rid = 2) as well as unrooted fermions with nt — 8 (n u = = 4). 

For the n t — 2 runs the gauge coupling was (3 = 4.8. The lattice spacing was a — 0.41 fm, measured from 
the string tension a using the value of \[a = 465 MeV [11]. For the nt = 4 case (i = 4.3 and a = 0.42 fm and 
for the nt = 8 case (i = 3.8 and a = 0.44 fm. In all three cases the bare quark mass was am q = 0.04 and the 
spatial extension of the lattice was L s = 6. In the two rooted case we used temporal lattice extensions of 
L t = 50, 100, 200, 300 while in the unrooted case only L t = 100, 200 was used. Table 1 contains the number 
of configurations for each setup. 

Using equation (45) the difference of the free energies aFj 3= i — aFf 3=0 were measured on each set of 
configurations. These are listed in Table 2. According to equation (11) the mass of the ground state in the 
h = 1 channel can be obtained using a linear extrapolation to T = 0. For comparison we measured the pion 
mass in all cases using the ordinary spectroscopic method, which will be denoted by m,r jS p. The measured 
free energy values, the linear fits and the comparisons to the spectroscopic pion masses can be seen in Figure 
3. 

As can be seen from Figure 3 the mass of the ground state in the ^3 = 1 sector agrees with the spectro- 
scopic pion mass within errorbars for both rooted and unrooted staggered fermions. 

8.2 Quenched case 

Equation (46) can be rewritten as 
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=i ;7r = lim 

L t — >oc 



L t V64 




3V 
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3V 



k=l 



(47) 
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Figure 3: The differences of the free energies of the isospin one and isospin zero sectors as a function of the 
temperature on dynamical staggered configurations with a spatial volume of 6 3 . The dashed lines show the 
linear fits to the data points. The T — > extrapolated values are compared to the spectroscopic pion masses. 



1/oT 


Number of configurations 


48 


2502 


96 


1852 


192 


731 


384 


412 



Table 3: Number of quenched configurations. 




The r.h.s. of equation (47) does not explicitly contain the number of staggered tastes. The quantity 

3V 2 \ 

E A * ) ( 48 ) 

fe=l / 

can be evaluated on quenched configurations as well. The question arises naturally: If one measures the 
pion mass on a quenched ensemble using regular staggered spectroscopy and evaluates the expression in (47) 
with the same fermion mass, will these be the same? 

To find this out we performed calculations on quenched configurations generated using the Wilson pla- 
quette gauge action. The spatial extension of the lattice was L s = 6, the gauge coupling was (3 = 5.6 
and the corresponding lattice spacing was a = 0.21 fm [11]. The time extension of the used lattices were 
L t = 48, 96, 192, 384 and for the measurements we used a bare quark mass of am q — 0.04. The number of 
configurations used are listed in Table 3. The results are summarized in Table 4 and the linear extrapolation 
is shown in Figure 4. The comparison shows that the result obtained from the free energies is consistent 
with the spectroscopic pion mass. 

The partition function contains all the information about the degrees of freedom present in the system, 
therefore, the free energy should be able to make a difference between dynamical and quenched configurations. 
The results, however, show that both types of ensembles yield a free energy that is consistent with particles 
of mass equal to the spectroscopic pion mass present in the system. Thus, one cannot tell this way whether 
a given set of configurations is from a dynamical or a quenched ensemble. 



9 Conclusions 

We have proposed a spectroscopic method based on purely thermodynamical considerations. The formulae 
obtained show the relation between the eigenvalues of the reduced staggered fermion matrix and the hadron 
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1/aT 


"of>,=i - aF /3=0 " 


48 


0.5393(6) 


96 


0.5389(2) 


192 


0.5385(3) 


384 


0.5389(4) 


— > oo 


0.5385(3) 




0.53874(3) 



Table 4: The "differences of the free energies" of the isospin one and isospin zero sectors, their T — » 
extrapolated value and the spectroscopic pion mass on quenched configurations. 



0.54 



0.539 - 



0.538 




Figure 4: The "differences of the free energies" of the isospin one and isospin zero sectors as a function of 
the temperature on quenched configurations with a spatial volume of 6 3 and bare quark mass am q = 0.04. 
The dashed line shows the linear fit to the data points. The T — > extrapolated value is compared to the 
spectroscopic pion mass. 
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spectrum. The method not only clarifies the findings of Ref. [1] in connection with the Goldstone pion mass, 
but also extends them. In principle, the method can be used to obtain the mass of the lightest particle 
in a given quark number sector. For example, in principle, the mass of the di-baryon could be obtained. 
However, it turns out that the application even to one-baryon states is computationally very demanding. 
Nevertheless, we successfully applied our method to the Goldstone pion. In the calculations presented the 
mass of the lowest state in the I3 = 1 sector is in agreement with the pion mass obtained using the ordinary 
spectroscopic method. This indicates that the method presented in Section 7 is a valid way of finding the 
pion mass. 
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